UWO-TH-09/11 



Causality of Holographic Hydrodynamics 

Alex Buchel 1,2 and Robert C. Myers 1 

1 Perimeter Institute for Theoretical Physics 
Waterloo, Ontario N2L 2Y5, Canada 

2 Department of Applied Mathematics 

University of Western Ontario 
London, Ontario N6A 5B7, Canada 

Abstract 

We study causality violation in holographic hydrodynamics in the gauge theory/string 
theory correspondence, focusing on Gauss-Bonnet gravity. The value of the Gauss- 
Bonnet coupling is related to the difference between the central charges of the dual 
conformal gauge theory. We show that, when this difference is sufficiently large, causal- 
ity is violated both in the second-order truncated theory of hydrodynamics, as well as 
in the exact theory. We find that the latter provides more stringent constraints, which 
match precisely those appearing in the CFT analysis of Hofman and Maldacena. 
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1 Introduction 

Hydrodynamics organizes the description of the macroscopic evolution of systems in 
local, but not global, equilibrium in terms of a derivative expansion. For concreteness, 
we will consider a four-dimensional relativistic fluid here. In the simplest situation 
(with no conserved charges), the dynamics of the hydrodynamic fluctuations in the 
fluid is simply governed by conservation of the stress-energy tensor T^ u , 

V^ = 0. (1.1) 

The stress-energy tensor includes both an equilibrium part (with local energy density 
e and pressure P) and a dissipative part Yi^ u , 

r^euY + PA^ + n^ where A^ u = + u^u v . (1.2) 

Above, is the local four- velocity of the fluid with u^u^ — — 1. Further, IP"^ = 0. 
In phenomenological hydrodynamics, the dissipative term W can be represented as 
an infinite series expansion in velocity gradients (and curvatures, for a fluid in a curved 
background), with the coefficients of the expansion commonly referred to as transport 
coefficients. The familiar example of the Navier-Stokes equations are obtained by 
truncating IP" at linear order in this expansion 

W u = nf (77,0 = -^""-(A^V^u, (1.3) 
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where 




Notice that at this order in the hydrodynamic approximation we need to introduce only 
two transport coefficients, namely the shear 77 and bulk ( viscosities. In the following 
discussion, we will be particularly interested in describing conformal fluids, in which 
case we must also impose the vanishing of the trace of the stress tensor. This restriction 
fixes C — 0, as well as P = e/3 in the equilibrium contribution. 

As noted above, hydrodynamics can be regarded as giving a systematic derivative 
expansion. Within this framework, it is then straightforward to extend IP^ to the 
next order to including terms of order V 2 w or (Vu) 2 . In general, this extension would 
require the introduction thirteen new transport coefficients [1]. However, if we again 
restrict our attention of conformal fluids, the second-order term Yl 2 only depends on 
five of these new transport coefficients [2]. While the interested reader can find the 
complete description in [2], we only illustrate the extension here by showing the first 
few new terms: 



The terms controlled by A2,3 involve the vorticity while k term is proportional to the 
spacetime curvature. Hence the terms explicitly given above are sufficient to describe 
the vorticity-free flow of a conformal fluid in a flat background spacetime. As noted 
in [2], the first term proportional to Tn essentially captures the second-order formalism 
of Miiller, Israel and Stewart (MIS) [3] while the subsequent nonlinear terms already 
represent an extension of their approach. However, this linear term is sufficient to ad- 
dress the question of causality within the hydrodynamic framework. It is well known 
that if the dissipative contribution is truncated as in (1.3), for any viscosity coefficients 
{77, (}, there are always linearized fluctuations for which the wave-front speed is super- 
luminal [4] . The primary motivation of MIS was then to eliminate this acausality in the 
hydrodynamic equations. Indeed the MIS term is sufficient to tame the superluminal 
propagation with an appropriate choice of the relaxation time Tn, as we demonstrate 
below [5]. However, we add that, as will become evident, the constraints on Tn emerge 
from the behaviour of modes outside the regime of validity of the second-order hy- 
drodynamic framework, i.e., from very short wavelength modes. Hence, one should 
keep in mind that these constraints do not signal any fundamental pathologies but 




(1.5) 



= —rj a 1 




3 



rather only indicate where a certain approximate mathematical framework describing 
the fluid becomes problematic. Nevertheless, a causal system of second-order hydrody- 
namic equations is still required in many situations, such as, numerical simulations [6] 
which implicitly extrapolate the hydrodynamic equations to the smallest numerical 
scales, even though the physics of interest is in the long wavelength regime. 

Linearized fluctuations of the second-order truncated hydrodynamics in conformal 
fluids were discussed in [2]: 

■ The dispersion relation of the shear channel fluctuations is given by (see eq. (3.27) 
of [2]) 



-ro 2 r n T-^ + fc 2 ^ = 0, (1.6) 

Z7T S 

where tt> = u/(2irT) and k = k/(2nT). Now the speed with which a wave- front 
propagates out from a discontinuity in any initial data is governed by [7] 



Re fro ) 

hm 



[shear] 



|fc|^oo k 

Hence causality in this channel imposes the restriction 

mT > - . 

s 



V = front /, ? x 

U [shear! • K 1 -') 



s r n T t shcar ] ' 



■ The dispersion relation of the sound channel fluctuations is given by (see eq. (3.20) 
of [2]) 

q m iw 2 1 m , o 4r? o ik 2 
-m 3 r n T - — + -mT wfc 2 + vok 2 + — = . 1.9 

2n 3 3s 07T 

Hence 

Re fro 



lim 



1 , 4?7 1 Jront 



[sound] v 3 3 S r n T ^ " 



|fc|^oo k 

From (1.10), causality in the sound channel imposes the following condition 

rnT >2^. (1.11) 

s 

From (1.7) and (1.10) above, we might note that both f ^°ar] an< ^ ^jsound] diverge as r n — > 
0. We may also see that the front velocity in the sound channel is always larger than 
that in the shear channel 1 and hence the former provides a more stringent constraint 
(1.11) on the transport coefficients of the second-order hydrodynamics. Again, we note 



1 In fact, this is a general result which extends to nonconformal fluids as well [8]. 
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that as should be evident from (1.7) and (1.10), these restrictions arise from pushing 
the second-order hydrodynamic framework beyond its natural regime of validity i.e., 
\k\ <C 1. We return to this point in greater detail in section 4. 

In principle, all the transport coefficients are determined by parameters of the un- 
derlying microscopic theory. In practice, such computations are prohibitively compli- 
cated as one has to derive the effective theory of hydrodynamics for a given microscopic 
system. The difficulties become even more insurmountable for strongly coupled plas- 
mas, as might be of interest at RHIC (or the LHC). However, the AdS/CFT correspon- 
dence of Maldacena [9, 10] provides a new framework in which transport coefficients 
are readily calculable at least for certain strongly coupled gauge theories [2, 11, 12]. 
Furthermore, it is a context where the discussion of conformal fluids becomes partic- 
ularly relevant. With reference to the constraints above, one has r]/s — 1/(Att) [13] 
and tjiT = (2 — log2)/(27r) [2,11] for N = 4 super- Yang-Mills or any strongly coupled 
four-dimensional gauge theory for which the holographic dual is described by Einstein 
gravity [14]. Hence a second-order hydrodynamic analysis of such holographic plasmas 
does not suffer from any problems with acausality. 2 

In this paper we extend this analysis using a particular effective model in the gauge 
theory/string theory correspondence. Specifically, we consider a holographic model 
with a Gauss-Bonnet (GB) gravity dual, 



1 = W P ' d X ^ 



12 A 

— + R + -^L 2 {R 2 - 4R, U R^ + R^R*""") 



1.12) 



The corresponding conformal gauge theory is distinguished by having two distinct cen- 
tral charges [16, 17] - see section 4 for details. The effect of such curvature-squared 
interactions on holographic hydrodynamics was examined in the context of string the- 
ory in [18], however, only within a perturbative framework. The GB gravity theory 
(1.12) is particularly well-behaved allowing the holographic analysis to be extended to 
finite values of the coupling A GB . In particular, the ratio of the shear viscosity to the 
entropy density is found to be [19] 



rj 1 

S 4lX 



1 -4A r 



b] ■ (1.13) 



Below we compute the relaxation time of the CFT plasma dual to GB gravity (1.12), 
Tn = 7T[(Agb). We find that the causality condition (1.11) then constrains A GB both 



2 A full analysis of causality in these holographic plasmas was also discussed in [15]. 
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from above and below. Note that in this case, the constraints are imposed to avoid 
fundamental inconsistencies in the theory. Further, these constraints on A GB would in 
turn lead to bounds on the viscosity in GB hydrodynamics. 

The analysis of the second-order hydrodynamics in GB gravity is interesting because 
causality violations were already used to produce an upper bound on the GB coupling 
in [20]. The analysis there also examined the propagation of signals through the dual 
gauge theory plasma but made no restriction to second-order hydrodynamics. Hence we 
turn to the study of causality violation in an exact analysis of the GB theory in section 
3. Following [21], the dispersion relation of physical fluctuations in a gauge theory 
plasma is identified with the dispersion relation of the quasinormal modes of a black 
hole in a dual gravitational description. Extending analysis of [20], we study dispersion 
relation of the GB BH quasinormal modes in the shear and the sound channels. As was 
done for the scalar channel in [20], we show that requiring that these modes are not 
superluminar, i.e., the phase velocity remains less than one in the infinite momentum 
limit, constraints A GB . We find that combined these constraints are more stringent than 
the causality constraints coming from the second-order truncated GB hydrodynamics. 



2 Causality of second-order Gauss-Bonnet hydrodynamics 

We are interested in determining when the second-order hydrodynamics dual to GB 
gravity satisfies the causality constraint (1.11). Since the ratio of shear viscosity to 
entropy density is already given in (1.13), it only remains to determine the relaxation 
time Tn for the dual plasma. The simplest approach to discover rn(A GB ) is to examine 
the dispersion relation of the sound quasinormal mode of a GB black hole. Here, the 
field theory considerations establish that [2]: 

ro = c s k - 2ni TT k 2 + ^flL ^ TrT _ I rT ^J k 3 + q(jA) . (2.1) 

With 

rr = l a,ld C ' = 7T (2 ' 2) 

this expression is simply the Taylor series solution of the dispersion relation (1.9). 

In the dual gravitational description the dispersion relation (2.1) is obtained by 
imposing the incoming wave boundary condition at the horizon and the Dirichlet con- 
dition at the boundary on the sound channel quasinormal mode wavefunction. The 
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technique is clearly explained in [21]. Here, we only present the salient steps in our 
analysis 3 . 

The planar black hole solution in GB gravity can be written as [22,23] 



r 2 



where 



and 



/M= W1-4A ,(1-^ (24) 

A 2 = \ (l + a/1 — 4A GB j . (2.5) 

The horizon is located at -u = 1 and asymptotic boundary is reached with u — > 0. 4 
Note that the normalization constant A is chosen so that A 2 f(u = 0) = 1. Hence the 
asymptotic behaviour of the metric shows that the AdS curvature scale is AL. The 
Hawking temperature, entropy density, and energy density of the black hole are 

Now the sound channel quasinormal mode satisfies the following equation 

Zfsound] ( U ) + ^sound ^[sound] ( U ) + ^sound ^[sound] (u) — , (2.7) 

where the coefficients C^ und are presented in Appendix A. In the hydrodynamic limit, 
k, to — > with ^ kept fixed, the incoming wave boundary condition at the horizon 
implies 

Z[sound] = (1 - u 2 y % ~% ^z (u] to, fe) + ik zi(u; to, fe) + k 2 z 2 (u; to, k) + C(fc 3 )^ , (2.8) 



with 



Zi(u; \i to, \i k) = Zi(u; to, k) , for any fj, , i = 0, 1, 2 



o 



lim to, fe) = 5 

U— >1 



(2.9) 



3 Further computational details are available from the authors upon request. 

4 A more conventional radial coordinate [19,22] would be given by r 2 = r\ju. Implicitly we have 
also chosen the branch of well-behaved solutions and we are restricting our considerations to A GB < 1/4 
- see [22,23] for details. 
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In (2.8) we kept terms in the hydrodynamic expansion to the order necessary to identify 
(2.1). The sound wave dispersion relation (2.1) is then obtained by imposing the 
Dirichlet condition on Z[ SOU nd] at the boundary 

limZ sound = 0. (2.10) 

To leading order in the hydrodynamic approximation we find 

_ k 2 (x + 4A GB - 1)A 2 - QX GB W 2 x 
Z ° 2X GB (2A 2 k 2 - 3n 2 )x ' 1 ' j 

where we used a more convenient radial coordinate 

(1-4A gb + 4A gb m 2 ) 1/2 . (2.12) 



x 



Imposing the Dirichlet condition (2.10) at this order recovers the conformal sound 
speed: c s = l/y/3. 

To order O(k) in (2.8), we find 

Zl =£- (\ GB (2A 2 k 2 - 3ra 2 ) x) x ( 2{A 2 k 2 {x + 4A GB - 1) - 6ra 2 A GB x) In 



8k x 

+ (1 - x)A 2 k 2 (x 2 -6x + 36X GB x + 12A GB - 3) + 6tv 2 X GB x(x + 3)(x - 1) 

(2.13) 

Imposing the Dirichlet condition (2.10) at order 0(k 2 ) identifies 

1 = J_[i_4A GB ], (2.14) 

S 47T L J 

in precise agreement with the result (1.13) which was originally determined with a 
Kubo formula computation in [19]. 

Unfortunately, we were not able to evaluate Z2 (and as a result rn) analytically. 
Thus, we had to resort to numerical analysis. For the question of causality, we are not 
interested here in rn per se, but rather in the relation (1.11). Hence figure 1 presents 
the difference (tuT — 2^) as a function of A GB . We find that unless A GB G [A m j„, A mox ], 
where X m i n = —0.711(2) and X max = 0.113(0), causality of the second-order truncated 
hydrodynamics of the GB plasma is violated. In figure 2, we also present the front 
velocities in the shear (1.7) and sound (1.10) channels. Note that we find that the 
relaxation time vanishes for A GB = 0.165(5), which causes both of the front velocities to 
diverge at this point in the figure. The graph also demonstrates that v^^r] < v [sound\i 
as noted above. 
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Figure 1: (Colour online) Causality of the second-order Gauss-Bonnet hydrodynamics 
is violated once TnT < 2 2 . Thus, A GB G [A min , X max ), where X min = —0.711(2) and 
X max = 0.113(0). 

3 Causality of full Gauss-Bonnet theory 

In the previous section we showed that depending on the value of a the GB coupling, 
A GB , or equivalently on the values of the microscopic parameters in the dual CFT, 
causality of the second-order truncated hydrodynamics of the dual theory can be lost. 
In this section, we wish to compare those results to causality violations found with the 
analysis of [20]. While the latter also looks for the appearance of superluminal signals 
propagating in the dual plasma, it makes no reference to truncating the derivative 
expansion in the hydrodynamic framework. We will find that the constraints on A GB 
arising from the exact analysis are more restrictive than those found in the truncated 
second-order hydrodynamic analysis in the previous section. 

Dispersion relation of the linearized fluctuations in plasma is identified with the 
dispersion relation of the quasinormal modes of a black hole in the dual gravitational 
description. There are three types of quasinormal modes in gravitational geometries 
with translationary invariant horizons [21,24]: 

■ a scalar channel (helicity-two graviton polarizations); 

■ a shear channel (helicity-one graviton polarizations); 
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front 




Figure 2: (Colour online) Front velocity for the shear (red) and sound (blue) channels 
for the second-order hydrodynamics, as given in (1.7) and (1.10). The dashed vertical 
lines indicate \ m in and X max , where vfj°™* d] reaches one. 



■ a sound channel (helicity-zero graviton polarizations). 

While the shear and sound channels correspond to those considered in the previous 
discussion of second-order hydrodynamics, the scalar channel was not mentioned there 
because it contains no modes whose frequency vanishes as k — > 0. Of course, this is in 
agreement with the standard hydrodynamic analysis [4]. However, the scalar channel 
quasinormal modes of the GB black holes in the limit k — > oo were studied in detail 
in [20]. It was found there that requiring 



lim — 

fc-^oo k 



<i, (3.1; 

[scalar] 



constraints A GB as follows 



A GB < KT r = ^ • (3-2) 

Note that X s ^ lar < X ma x found in the context of the second-order truncated hydrody- 
namics. In the remainder of this section we extend analysis of [19,20] to the shear and 
the sound channel quasinormal modes. 
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3.1 Causality in the shear channel 

It is straightforward to derive the shear channel quasinormal equation for the GB black 
holes: 

Zfshcar] («) + C^ £ar Z[ shcar] («) + Cfl ar Z [shcar] (u) = , (3.3) 

where the coefficients C^ ear are presented in Appendix B. Following [20], we now caste 
(3.3) into the form of the Schrodinger equation. Towards this end, we introduce a new 
radial coordinate y 



dy 



du ■u 1 / 2 /(m) ' 

and rescale the radial profile as 



(3.4) 



Z [shear] — ^ 1p [shear] , (3-5) 

with 

^ In B = (^4 M (2A GB / - 1) V(2A GB / - l) 2 + -4 2 /(4A GB - 1))) 
x ^ 2 (1 - 4A GB )(12/ 3 A 2 B - 16/ 2 A GB - 4 + 7f) - (20/ 2 A 2 B - 20A GB / (3.6) 
+ 8A GB + 3)(2A GB /-1)V^) , 



where a = xo/k. The quasinormal equation (3.3) can then be rewritten as 

— h 2 dy ^[shear] + ^[shear] ^[shear] = O? ?/>[ s hear] , ^ = T ) 

« (3.7) 

where <7 [shcar] = U^ sheax] + f? C/ [shear] . 

The first part of the effective potential has the simple form when expressed in terms 
of u 



"[shear] («) " _ ^2 



(1-4A GB )(1- y/l-4A GB (l-^)) 
(1 - 4A GB (1 - u 2 )) (1 - VI - 4A GB ) ' 

while the expression for U] hear is too long to be presented here, but we note that the 
latter is a function only of u, A GB and a. What is important is that in the limit k — > oo 
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Figure 3: (Colour online) Typical behaviour of U°, the leading contribution to the 
Schrodinger potential, for both the shear (blue) and sound (red) channels. The solid 
and dashed curves show the behaviour for large and small |A GB |, respectively. (Our 
representative values here are: A GB = —1.5 and —0.15.) 

(or H — > 0), everywhere except in the tiny region y > — ^ the dominant contribution 
to U shear comes from U® hear . Thus in this limit we simply replace 

o i fo y<0, 

Figure 3 illustrates the general behaviour of the leading potential (3.8). For any 
values of A GB , we have in the asymptotic region, U? ghBBX , (u = 0) = 1 while at the horizon, 
^[shear](' u = -0 = 0- ^ow f° r sma U values of |A GB |, ^hear] * s a monotonically decreasing 
function between these two points. However, for larger negative values of A GB , the 
potential develops a (single) maximum at intermediate value of u: 

4(V1 - 4A OB - 1) 2A GB v ' 

As might be inferred from u max above, the critical coupling for the appearance of 
this maximum is A GB = —3/4. At this stage, the analysis is identical to that for 
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the scalar channel studied in [20]. Once the effective potential in the Schrodinger 
problem (3.7) develops this new maximum, there always exist quasinormal modes with 
Re (a 2 ) ~ U^ ax > 1. This implies then that in the limit of infinite fc, Re(tt>)/fc > 1 
for these modes and hence they lead to a violation of causality. Hence requiring the 
excitations in the shear channel to be well behaved imposes the constraint: 

a gb > K h r = -\- (3.H) 

3.2 Causality in the sound channel 

The quasinormal equation for the sound channel is given in (2.7). Following [20] (as 
reviewed in the previous section), we arrive at the corresponding Schrodinger problem 
in the sound channel 

- h 2 C^^Vund] + C/[ sound ] ^[sound] = O? ^[sound] • (3-12) 

Once again, in the limit k — > oo, potential £/[ soun d] is given by 

Tj J ^[sound] 2/ <0 > /oioN 

^ [sound] = < (3.13) 

I +00 y — . 

where 

TT o (l-8A GB -4/A GB (A GB /-l))A 2 / 

^[sound] - (2A GB / - 1) 2 (3 - 14) 

(1 - 4A GB (1 + u 2 )) (1 - y/1 - 4A GB (1 - it 2 ) ) 
(1 - 4A GB (1 - « 2 )) (1 - VI - 4A GB ) 

The general behaviour of this potential is the same as described in the previous section, 
as can be seen in figure 3. In particular, to avoid the appearance of an intermediate 
maximum in the potential (3.14) and the corresponding causality- violating quasinormal 
modes, we must impose the constraint 

A GB >A G °r d = -^. (3.15) 

Note that A G ° B W > A G g ear found above and also \ s °£ nd > \ min found in the context of 
the second-order truncated hydrodynamics. 
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4 Conclusion 



In this paper, we have shown how causality of the near-equilibrium phenomena in the 
quantum field theory can be used to distinguish "healthy" models from the "sick" ones. 
To illustrate the point, we studied the fluctuations in the CFT plasmas, holographically 
dual to Gauss-Bonnet gravity (1.12). Our analysis found two sets of constraints on the 
GB coupling: 



It is clear that the exact analysis of section 3 produced more stringent restrictions than 
the analysis of the truncated second-order hydrodynamic equations in section 2. While 
both of these approaches are examining the behaviour of gravitational fluctuations in 
the GB black hole background (2.3), the relevant quasinormal modes are very different 
in the two cases. The second-order hydrodynamics is focused entirely on the behaviour 
of the sound mode, i.e., the lowest quasinormal mode in the sound channel. In contrast, 
the potential causality violation by highly excited quasinormal modes in the scalar 
channel set the upper bound in (4.2) while the lower bound arises from a similar set of 
quasinormal modes in the sound channel. 

We must emphasize that the status of these constraints differs at a very basic 
level. Theories outside of the bounds given in (4.2) are fundamentally pathological. In 
contrast, the constraints (4.1) simply indicate where a certain approximate description 
of the fluid becomes problematic. As such, it is somewhat remarkable then that the 
bounds coming from these two very different approaches seem to be fairly close to 
each other. It is also satisfying that the fundamental constraints (4.2) are the most 
restrictive so that the truncated second-order hydrodynamics will be stable in any of 
the cases where the underlying theory is physically sound at a fundamental level. A 
priori, this does not seem to be required by any basic principles. 

We consider the second-order hydrodynamics and the behaviour of the lowest sound 
quasinormal mode in more detail in figures 4 and 5, which show results for A GB = —2.5 
- note that the latter is outside the "healthy" ranges in both (4.1) and (4.2). Using 
the truncated second-order equations only yields physically reliable results for k <C 
1, where the dispersion relation can be Taylor-expanded as in (2.1). These Taylor 
expansions for the phase velocity and the width, i.e., Re(tu)/fe and Im(tt>), keeping 
only the 0{k 2 ) terms are illustrated with the green curves in the two figures. On the 



second-order hydrodynamics: —0.711 < A GB < 0.113 , 
exact analysis: — ^ < A GB < . 



(4.1) 
(4.2) 
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Figure 4: (Colour online) The phase velocity Re(tn)/fc for A GB = —2.5: The blue 
line shows the behaviour of the lowest quasinormal mode calculated numerically. The 
red curve corresponds the second-order hydrodynamic approximation (1.9). The green 
curve corresponds the next order Taylor expansion (2.1) arising from the second-order 
dispersion relation. 

other hand, the causality analysis of section 2 treats the dispersion relation (1.9) as 
exact and the results are shown with the red curves. In this case as k — > oo, the phase 
velocity rises to Re(tt>)/fe = v^J^ ~ 1.126 and the width also reaches a finite value 
asymptotically 



As is evident in the figures, these "exact" dispersion relations only match the Taylor 
expansion (2.1) for small values of k. 5 

To emphasize the limitations of treating (1.9) as an "exact" dispersion relation, 
figures 4 and 5 also show numerical results for the behavior of the sound mode, i.e., 

5 In fact, it may seem that, in figures 4 and 5, the two sets of curves begin to separate at surprisingly 
small wave numbers (i.e., k ~ 0.05 to 0.10). This observation can be explained by examining the 
Taylor series to higher orders. We find roughly an expansion in [(1 — 4A GB ) 2 fc 2 ] rl for large Aqb and so 
with Aqb = —2.5, the Taylor series should only be expected to match the "exact" result for k <C 1/11. 



lim 

|fc|— >oo 




[sound] 




-0.047. 



(4.3) 
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Figure 5: (Colour online) The decay width Im(lt)) for A GB = —2.5: The blue line 
shows the behaviour of the lowest quasinormal mode calculated numerically. The red 
curve corresponds the second-order hydrodynamic approximation (1.9). The green 
curve corresponds the next order Taylor expansion (2.1) arising from the second-order 
dispersion relation. 

the lowest quasinormal mode with the blue curves. As expected, all of the different 
curves agree at small k but not at large k. While not unexpected, we wish, in particular, 
to point out that the striking differences between the actual behaviour of the sound 
mode and the second-order hydrodynamic dispersion relation (1.9). Figure 5 shows 
that the actual width decays rapidly to zero in contrast to the finite asymptotic limit, 
given in (4.3) above, for the second-order dispersion relation. Similarly in figure 4, 
the actual phase velocity rises beyond i^^] — 1.126, the asymptotic limit found for 
second-order hydrodynamics, but then appears to decay back towards one as k — » oo. 
While the numerical results shown are already becoming less reliable for k > l, 6 it 
seems that the limit k — > oo should produce a front velocity which respects causality. 7 
Despite the fact that the results of (1.9) may have little resemblance to the physical 
behaviour of the sound mode at large k, it remains important that the truncated 

6 These numerical difficulties are correlated to the dramatic decrease in |Im(ro)|. 
7 Of course, the higher quasinormal modes are expected to violate causality with v max ~ 1.415, 
using the analysis of the effective Schrodinger potential (3.14). 
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second-order hydrodynamic equations present a robust mathematical framework in 
certain situations. For example, numerical simulations of the strongly coupled quark- 
gluon plasma [6] implicitly extrapolate the hydrodynamic equations to the smallest 
numerical scales, even though the physics of interest is in the hydrodynamic regime. 
In this situation, having a system of hyperbolic equations which is causal is of course 
essential. 

Our holographic construction with Gauss-Bonnet gravity provides a simple toy 
model in which the dual CFT is completely specified by two central charges c and a. 
The exact relation of these CFT parameters to the gravitational couplings in the action 
(1.12) is given by [16,25] 

2 r3 



TT 2 L 



GB ; 

2 r3 



2 3 /2 £3 p 



F (1 + V1-4AOB) 3 / 2 (3y/l-4\ GB - 2) , (4.4) 



and hence 



- - 1 . (4.5) 

V1-4A GB J 

As A GB alone fixes this last combination, we may re-express our causality constraints 
(4.2) in terms of the central charges, which yields the elegant result: 

-I < ^ < I . (4.6) 
2 - c ~ 2 K J 

Here we have focused on the fundamental constraint (4.2) rather than the second-order 
hydrodynamic constraint (4.1) and have found that if the difference in the central 
charges grows too large, some linearized fluctuations in the model propagate faster 
than the speed of light. 

These constraints are also intriguing in comparison to the analysis of four-dimensional 
CFT's presented by Hofman and Maldacena [26]. They consider "experiments" in 
which the energy flux was measured in various directions at null infinity after a lo- 
cal disturbance was created in the stress energy. This energy flux was found to be 
controlled by the three-point function of the stress tensor in the CFT. Further it was 
shown that the parameters fixing this three-point function must be constrained in order 
that the energy flux was always positive. The holographic dual describing this CFT 
"experiment" would involve both curvature-squared and curvature-cubed interactions. 
However, if the CFT was restricted so as to eliminate the curvature-cubed interaction 
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in the gravitational dual, 8 then the positive energy constraint reduced to a constraint 
on the central charges, which in fact precisely matches that given in (4.6). The precise 
agreement of the upper bound was already noted in [26]. 

One can go further in that measurements in the above "experiment" can be or- 
ganized into three different channels, just as for the graviton fluctuations in section 
3. With this classification, the upper bound in (4.6) comes from the appearance of 
negative energy flux in the scalar channel in [26] while it arises from causality violation 
in the same set of fluctuations in the holographic calculations [20] . Similarly, the lower 
bound is set to avoid problematic behaviour in the sound channel in both approaches. 
Then we may also note that while it did not set a fundamental constraint, causal- 
ity violations also appear in the shear channel at the critical value given in (3.11): 
^GB £ar — — f ■ This result then translates to a critical value (c — a)/c\ s h ea r — — 1, which 
again precisely matches that for the appearance of negative energy fluxes in the shear 
channel. Hence, at least within this holographic model, we have drawn a precise cor- 
relation between the appearance of negative energy fluxes and of superluminal signals 
in various channels. 

We close with a few comments about other potential instabilities in this holographic 
model. It was observed in [19] that a new instability arises in the dual plasma at 
A GB = —1/8. At this point, the effective Schrodinger potential develops a small well 
where U° < just in front of the horizon (i.e., near u — 1). For sufficiently large k 
(and |A GB |), this well will support unstable quasinormal modes, as described in [27]. 9 
Examining (3.14) reveals similar behaviour and hence instabilities for A GB > 1/8 in 
the sound channel, however, this problem only appears outside of the range (4.2) 
allowed by causality. Formally, the effective potential (3.8) in the shear channel shows 
a similar behaviour for A GB > 1/4 but this is again of no consequence since, as noted 
before, our entire analysis is only valid in the regime A GB < 1/4. Hence the instability 
in the scalar channel is the only one that appears within the physical regime (4.2). 
This instability does not correspond to a fundamental pathology with the theory but 
rather indicates that the uniform plasma becomes unstable with respect to certain non- 
uniform perturbations. It would of course be interesting to follow the full nonlinear 
effect of these instabilities. On the gravitational side, this instability seems similar in 
certain respects to the Gregory-Laflamme instability for black strings [29]. 

8 This restriction would be realized in any supersymmctric CFT. 

9 It was later observed that this instability seems to be increased by a chemical potential [28] . 
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A Coefficients of (2.7) 



(i) 



fu (fA 2 k 2 + 2A 2 k 2 - 3ro 2 ) - 2f 2 u (fA 2 k 2 + lOA 2 k 2 - 12ro 2 )A, 



sound 



■GB 




'GB 



(A.l) 
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Cfl nd = (/V [-(fA 2 k 2 + 2A 2 k 2 - 3m 2 ) + 2 f(12A 2 k 2 + 2 fA 2 k 2 - 15ra 2 ) A c 

- 24/ 2 (-5ro 2 + 4^ 2 fc 2 )A 2 B - 16/ 3 (15ro 2 - 10^ 2 fc 2 + fA 2 k 2 )X 3 GB 

+16/ 4 (15m 2 - 6^ 2 fc 2 + fA 2 k 2 )X 4 GB - 96/ 5 ra 2 A GB ] ^ x ^3rv 4 u - AfA 2 k 2 

+ f 2 A 4 k 4 u + 8f 2 A 2 k 2 - AfA 2 k 2 - Atv 2 fA 2 k 2 u - 2n> 2 A 2 k 2 u + 2A 4 k 4 fu 
+ 2f(5fw 2 uA 2 k 2 - 8fA 2 k 2 - 15tv 4 u - 28fA 2 k 2 - 8A 4 k 4 u + 8A 2 k 2 
+ 2f 2 A 4 k A u + 24rv 2 A 2 k 2 u - 8fA 4 k% + 28/l4 2 fc 2 )A GB - 4/ 2 (-30ro 4 w 

- 12/ro V4 2 fc 2 - 32A 4 k 4 u - Af 2 A 2 k 2 + 2f 2 A 4 k 4 u + 12fA 4 k 4 u + A0fA 2 k 2 
+ 60m 2 A 2 k 2 u - 2AA 2 k 2 - 11 fA 2 k 2 )\ 2 GB - 8f{2fA 4 k 4 u - 56ra 2 A 2 k 2 u 

+ 3fU 2 fc 2 - 2AA 4 k 4 fu + 22f 2 A 2 k 2 + 26rv 2 fA 2 k 2 u - 24fA 2 k 2 + 30ro 4 w 
+ 24A 4 k 4 u)X 3 GB + !Qf 4 (lQw 2 fA 2 k 2 u - 18w 2 A 2 k 2 u - QA 4 k 4 fu + Qf 2 A 2 k 2 



+ f 2 A 4 k 4 u + 15tt> 4 u)\ 4 GB - 96 fu ra 2 (m 2 + fA 2 k 2 )\ 



GB 



(A.2) 



where / is given by (2.4). 



B Coefficients of (3.3) 

C£L = (-fu (-m 2 + fA 2 k 2 ) + Af 2 u {fA 2 k 2 - 2tv 2 + k 2 A 2 )\ c 



-i 



- 4/ d u (fA 2 k 2 - 6m 2 + 4fc^4 2 )A GB + 16 fu (-2tv 2 + A 2 k 2 )\ GB + lQfu K> 2 \ 4 C 
x (^f 2 A 2 k 2 - 2m 2 + fro 2 - 2f(4fA 2 k 2 + fro 2 - 4tt> 2 ) A GB + 8/ 2 (-/m 2 - tv 2 
+ 2k 2 A 2 )\ 2 GB + 24/ 4 A GB tx> 2 - W 5 A 4 GB m 2 ) 

(B.l) 

C2L = ((2A GB / - l) 2 A) ' x (vo 2 - fA 2 k 2 + Af(k 2 A 2 - m 2 )A GB + 4/ 2 A 2 B * 2 ) 

(B.2) 

where / is given by (2.4). 
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